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Abstract 



> 

en 

tH- ■ We discuss two physical examples of the so-called pseudo-bosons, recently introduced in 

f->-. ■ connection with pseudo-hermitian quantum mechanics. In particular, we show that the 

so-called extended harmonic oscillator and the Swanson model satisfy all the assump- 
tions of the pseudo-bosonic framework introduced by the author. We also prove that the 
biorthogonal bases they produce are not Riesz bases. 



I Introduction 

In a series of recent papers, [I]-[l], a family of excitations generalizing bosons have been intro- 
duced by considering an extended version of the canonical commutation relation, [a,a^ = 1, 
which look like [a,b] = 1, with b ^ a^. We have shown that, under suitable assumptions, two 
families of biorthogonal bases of the Hilbert space Ti on which a and b act can be constructed 
and that these sets may, or may not, be Riesz bases, see Section II. In particular, this character- 
istic is interesting since, as we have discussed in [1], any Riesz basis by itself produces non trivial 
examples of pseudo-bosons. In [3] we have shown how to construct examples of pseudo-bosons 
by generalizing standard techniques used in supersymmetric quantum mechanics, while in [1] 
we have analyzed some mathematical aspect of our construction. Here we continue our analy- 
sis more on a physical side, discussing in details the pseudo-bosonic structure for the extended 
quantum harmonic oscillator (EQHO), [5], and for the Swanson hamiltonian, [5l[6l[7]. 

The paper is organized as follows: in order to keep the paper self-contained, in the next 
section we quickly review the general framework of pseudo-bosons. In Section III we construct 
the biorthogonal bases for the EQHO. In Section IV we discuss the Swanson hamiltonian and 
its pseudo-bosonic excitations. Section V contains our conclusions. 

II The general settings 

In this section we will review the general framework originally introduced in [1] and further 
developed in [2l[3lll|. 

Let "H be a given Hilbert space with scalar product (., .) and related norm ||.||. We introduce 
a pair of operators a and b acting on "H and satisfying the following commutation rules 

[a,b] = l. (2.1) 

Of course, this collapses to the OCR's if 6 = a"'". It is well known that a and/or b are unbounded 
operators, so they cannot be defined in all of H. Following [1] we consider the following 

Assumption 1. there exists a non- zero (po ^H such that a(pQ = 0, and (fo G D°°{b). 
Assumption 2. there exists a non- zero \&o G ^ such that 6^^o = 0, and \&o G -D°^(a^). 

Here D^{b) is the domain of all the powers of b: D°°{b) = r\k>oD{b''), where D{b'^) is 
the domain of b''. Analogously D°°[a^) is the domain of all the powers of a^. Under these 



assumptions we can introduce the following vectors of l-L: 

</.„ = 4^6>o and ^„ = 4^(atr^o,, n > 0. (2.2) 

yrv. \/n\ 

Let us now define the unbounded operators N := ha and 91 := A^"'' = a)W . It is possible to 
check that ipn belongs to the domain of iV, D{N), and \E'„ G -D(9T), for all n > 0. Moreover, 

Nipn = riifn, 9I^„ = n^„. (2.3) 

Under the above assumptions it is easy to check that (\E'„, ipm) = 6n,m {^o, fo) for all n,m > 
0, which, if we chose the normalization of \i/o and (fo such that (\i/o, v^o) = I5 becomes 

{^n, ^m) = Sn,m, Vn, m > 0. (2.4) 

This means that the sets T^ = {\l/n, n > 0} and J-^^ = {(fn, n > 0} are biorthogonal and, 
because of this, the vectors of each set are linearly independent. If we now call V^p and V^ 
respectively the linear span of J^^ and J-"*, and Ti^ and "H^ their closures, then 

oo oo 

/ = ^(*„, /)(/.„, yfen^, h = J2{Vn,h)^rr, WheU^. (2.5) 

n=0 n=0 

What is not in general ensured is that all these Hilbert spaces do coincide, i.e. that "H^ = 7/^ = 
"H. Indeed, we can only state that "H^ C "H and Ti^ C "H. However, motivated by the examples 
already discussed in the literature and by the results in Sections III and IV, we consider 

Assumption 3. The above Hilbert spaces all coincide: T-dp = l-iq, = "H, 

which was introduced in [1]. This means, in particular, that both J-",^ and J^<!j are bases of "H. 
Let us now introduce the operators S^ and S*^ via their action respectively on J> and J^^: 

S^'^n = fn, S^ipn = ^n, (2.6) 

for all n > 0, which also implies that \I'„ = {S^iJ S^)'^n and c/?„ = {S^S<i,)(pn, for all n > 0. 
Hence 

Sq, Sip = Sip Sq, = 1. =^ Sqj = S^ . (2-7) 

In other words, both Sqj and S^p are invertible and one is the inverse of the other. Furthermore, 
we can also check that they are both well defined and symmetric, [T]. Moreover, if (\i/o, (po) = 1, 
it is possible to write these operators in a bra-ket language as 

oo oo 

S^ = J2\'^^><'^^\^ 5g, = ^|^n><^n|- (2.8) 

n=0 n=0 



These expressions are only formal, at this stage, since the series could be not converging and the 
operators S^ and S*^ could result unbounded. This aspect was widely discussed in P, where the 
role of Riesz bases in relation with the boundedness of S^ and Sq, has been discussed in detail. 
For completeness' sake we recall that, given an orthonormal (o.n.) basis Q = {(?„, ra > 0} of 
"H, and a bounded operator X with bounded inverse, Q^^^ := {gii := Xgn, n > 0} is a Riesz 
basis of "H, [8]. In particular, we have shown that S^p and S^ are bounded if and only if J^q, and 
J^ip are Riesz bases. In this case we will call our excitations regular pseudo-bosons, otherwise we 
will just call them pseudo-bosons. We will come back on this aspect in the following sections. 

We end this short review by recalling that these pseudo-bosons give rise to an interesting 
intertwining relation among non self-adjoint operators, see [I]-[l] and references therein. In 
particular it is easy to check that 

S^N = mS^ and N S^ = S^m. (2.9) 

This is related to the fact that the spectra of A^ and 91 coincide, and that their eigenvectors are 
related by the operators S^ and S^i,, in agreement with the literature on intertwining operators, 

[a [10]. 

Ill The extended quantum harmonic oscillator 

The hamiltonian of this model, introduced in [5], is the non self-adjoint operator Hj^ = 
I (p^ + x^) -|- iV2p, where /3 is a positive parameter and [x,p] = i. This hamiltonian is not 
PT-symmetric but satisfies VHis = HIV, where V and T are the parity and the time-reversal 
operators. In [3] the right and left eigenvectors, as well as the spectrum of if^, were found. 
Here we will see how this model can be discussed in the framework of Section II. In particular 
we will show that these right and left eigenvectors are nothing but the eigenvectors of -ff/? and 
Hi, which are automatically biorthogonal because of (12. 4p . We will also prove that the two 
resulting sets are both complete in £^(M), but they are not Riesz bases. 

Introducing the standard bosonic operators a = 4^ (x + ^), a) = -j={x — ^), [a, a)] = \, 
and the number operator N = a^a, we can write Hp = j3N -\- [a — a^) -\- ^\ which, introducing 
now the operators 

^/3 = a-^, Bp = a^ + ^, (3.1) 

can be written as 

Hp = l3{BpAp + ^p\), (3.2) 



where 7^3 = ^^- It is clear that, for all /3 > 0, AJ^ 7^ 5^ and that [Afs, Bfj] = 1. Hence we have 
to do with pseudo-bosonic operators. This, of course, does not mean that the Assumptions of 
the previous section are necessarily satisfied; in [H [3] we have discussed examples of operators 
satisfying this commutation rule and for which some of the assumptions (or even all!) do not 
hold. However, this is not the case here: we will show that the Assumptions 1, 2 and 3 are 
indeed satisfied. 

First of all, to check Assumption 1, we need to find a non zero vector c/?q G "H such that 
Ajs^Q = and c/?q G D^{Bp). This is quite easy since condition ApLpQ = can be rewritten 
as aifQ = ^ v9q . Hence, v^q is a standard coherent state with parameter -1: 

00 o—k 

^f = [/(r^Vo = e-V^^^5^^V^., (3.3) 

fe=0 

where ipo is the vacuum of a, a^po = 0, and U{(3~^) = e?^° '""' is the unitary (displacement) 
operator appearing in the theory of coherent states, [H]. Moreover, the vectors ipk = ^^y=- (fo, 
/c > 0, in (13. 3 P belong to the usual ortonormal (o.n.) basis J-" for an harmonic oscillator with 
hamiltonian h = lo [a'^a + |). Of course, for all possible choices of /3, we have Hv^q || = Hv^oll = 
1. In order to check whether c/?q belongs to D^^Bp) we repeat similar calculations as in |3]: 
since U~^{p-^)BpU{P'^) = 5^/2, and since B'^^if^^^ = U{(3-^) {U^\/3~^)B,3U{(3~^))'' ifo, it is 
easy to check that ||i?oV9o || < kle^^^, k > 0. Hence ^Pq belongs to the domain of all the 
powers of B^, as required. This is crucial, since it implies that 

V.f = 4=S,>f, (3.4) 



is well defined for all n > 0. 



To check now the validity of Assumption 2 we first look for the solution of BI'^q = 0. This 
can be rewritten as a\l/Q = — -^ \l/o , so that the solution is Wq ' = ^q = f/(— /3~^)(y5o = 
U~^{l3~^)(pQ. Using now the equality U{l3~^)A'^U~^{(3~^) = A'a,2 we can also check that 
WiAlf'^lf^W < k\e^/P, k>0. Hence 



^f = ^=iAlr¥,'\ (3.5) 



is also well defined for all n > 0. General reasons discussed in [T] show that, calling Nf^ = BpAp 

iV/3^i^^=n^f, %^(f)=nvl;(f), (3.6) 



and % = A^j = A^pB\, since 



these vectors above are biorthogonal and, since (v^o ' ^o / ~ ^ "^/P ^ the following holds: 

(^i^^^L^^> = ^n,™e-/^^ (3.7) 

Remark:— We could remove the factor e"^/^ by changing the normalization of Lp^ and 
\E'q . We prefer to keep this normalization since it is standard for coherent states. 

Remark:— Since in the ^-representation the coherent states of a = 4= (x + ^) with eigen- 
values ±4 look like -^e"^^^^'^'^'^' , it is possible to deduce the following expressions for the 
vectors introduced so far, which we now write introducing explicitly the dependence on x: 



and 



^mu) = — 1^ (x-^- y^] e-^(-+^/«^ 



7rV4v/2"r2! y dx (3 
Not surprisingly, these functions coincide, except for some normalization constants, with the 
right and left eigenvectors in |5]. 

Let us now define the following sets of vectors: J-'^ = {fn, n > 0} and J-"^ = {\E'n , n > 
0}, their linear span V^ and T)^ , and the Hilbert spaces "H^ and T-L^ obtained taking their 
closures. 

Concerning Assumption 3, we will now prove that J-'!f is complete in "H. To show this, we 
first observe that U-^{(3~^)a^U{l3^^) = Bp, so that 



Hence, we can check that a generic vector / G "H is orthogonal to all the Lpn s if the following 
scalar products are all zero: (U{—(3^^)f, ( a^ + 1 1 1 (/?o ) = for all n > 0. This immediately 
implies that (U{—(3~^)f^ {a})'^LpQ^ = for all n > 0, and therefore that U{—(3~^)f = 0, since J-' 
is an o.n. basis of "H. Hence / = 0. 

Similar techniques also show that J-"^ is complete in "H. We conclude that Assumption 3 
holds true. 

To check whether J-'!f and J-'^ are Riesz bases or not it is convenient to introduce the 
following self-adjoint, unbounded and invertible operator: Vg = e'-""'"" ^^^. This is useful because 
it turns out that 

Ap = VpaV^\ Bp = Vpa%\ (3.8) 
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It should be stressed that these equahties, as well as many of those which will appear in the rest 
of this section, are only defined on the dense set J-", since they involve unbounded operators. 

Formula fl3.8p implies that iJ^ can be related to a self adjoint operator hp = (3{a^a + 7/3I) 
as Hj3 = Vphj^Vp^ or, equivalently as 

HpVp = Vphp, (3.9) 

which shows that Vg is an intertwining operator (10) relating hp and Hp. Moreover, taking 
the adjoint of (13. 9p . we get VpWp = hpVp, so that Vp is also an 10 between hp and H\. This 
has well known consequences on the spectra of the three operators /1/3, Hp and H\ and on their 
eigenstates, [9]. 

In particular we have hpipk = e^ Vfc; where e^ = /3(fc + 7^3), V/c > 0. Hence, calling 
$^ := Vg (y9fc, which is clearly different from zero for all k since ker(Va) = {0}, we have, using 
the intertwining equality (13. 9p . 

Moreover, because of ( 13. 6p . since H^ip], = l3{Np + 7/3)v2fc = ef \r^\ and assuming that the 
eigenvalues ef are all non degenerate, it turns out that ipf = au $^ for all A; > 0, where ak 
are simply complex constants. As a matter of fact we can further check that all these constants 
coincide: a^ = e~^^^ , /c > 0, so that, in conclusion, 



^f)=e-V/5 Vp^u. (3.10) 



for all A; > 0. Similar arguments can be repeated to analyze the eigensystem of Hi: since 



Hi = V^^hpVp, if we put Xf^ = V^Vfe for all A; > 0, we get if J Xf ^ = ef Tf\ Moreover, 
using (13. 6p . H^ Wj^' = /3(Dl/3 + 7^)\l/[.^ = e[. ^k ■ Therefore, because of the non degeneracy of 
e[. , we again deduce a proportionality between T[.^ and \I/[. , which produces the equality 

*f = e-^/'^V.-V., (3.11) 

for all k > 0. This equation, together with (I3.10p . also implies that \1/^ = VT^y^jf , \fk > 0. 
Therefore, recalling (12. 6p . we recover the explicit expressions for the operators S^ and S^ . 
S^ = Vp"^ and, consequently, Slf = V^. This is in agreement with the following computations, 
extending formulas (12. 8p to the present situation where (y^Q , \1/q ) = e'"^/^ ^ 1: 



..00 / 00 \ 



V^O ' ^0 / n=0 \n=0 



,</?o ) n=0 \n=0 / 

where we have used the closure relation for T: XlfcLo Iv^™) (v^nl ~ ^- ^^ ^^^° have 



Wq ,(y9o ) „=o \n=0 / 

Going back to the nature of the sets T)p and T^\ since their vectors are obtained by the o.n. 
basis T via the action of the two unbounded operators V^ and Va^-, they are not Riesz bases. 
Hence, Assumption 4 in [1] is not satisfied: we have pseudo-bosons which are not regular. 

IV The Swanson hamiltonian 

The starting point is the following non self-adjoint hamiltonian, [5]: 

Ee = \ if ^^')-\ tan(2^) (p^ - x") , 

where 6' is a real parameter taking value in (— f , f ) \ {0} =: /. It is clear that Hq = H^g ^ Hq, 
for all 6 E I. As usual, [a;,p] = il. Introducing the annihilation and creation operators a and 
a^ we write 

He = N+'^ tan(2^) (a' + {a^f) + ^ 1, 

where N = a^a. This hamiltonian can be rewritten by introducing the operators 

Ag = cos{9) a + i sin(^) a\ 
Bg = cos{9) a^ + 2 sin(^) a, 

as 

Hg = Ug(BgAg + ^l\, (4.2) 

where ujg = ^^^2g) is well defined since cos(26') ^ for all 9 E I. It is clear that A'g ^ Bg and 
that [Ag, Bg] = 1. Hence we can try to see if it is possible to construct pseudo-bosons out of 
this hamiltonian (i.e. if Assumptions 1, 2 and 3 hold), or, even more, if regular pseudo-bosons 
arise (i.e. if also Assumption 4 of [1] is satisfied). It will be convenient to rewrite (14. ip by using 
the coordinate expressions for a and a^: 

^^ = 71 (^^'^ + '"' ^) ' u 3) 



We are now ready to check the vahdity of Assumptions 1 and 2. To begin with, we consider an 
abstract point of view: (^q satisfies AgLfQ = if and only if ay^Q = —i tan(6') a'^ (^q . Expand- 
ing 9?Q in both sides of this equahty in terms of the o.n. basis J-' = {yj^, n > 0} introduced in 
the previous section, ip^ = Xl^o "^nVn, and taking the scalar product of the resulting equation 
with ifj, j = 0, 1, 2, . . ., we deduce that C2n+i = and C2n = (— ztan(6'))"' y (ln)\\" ''O' '^^ — ^' 
Here (2n)!! = 2n-2(n-l) ■2(n-2) ■ ■ ■4-2 and (2n + l)!! = (2n + 1) ■ (2n- 1) ■ (2n-3) ■ ■ -S-S- 1. 
Therefore we get 



oo 



(2n-l)!! 



This vector belongs to "H for all 6 E I. Indeed we have ||(/9q p = |cop Yl'^=o (tan(6'))^ .2^^,, , 
which is always convergent for all 6' G /. 

Similar computations can be repeated to find the vector \I^q satisfying Bq^^q =0. In this 
case we find 

vl/(^) = rfo E (^ tan(^))« \ r'^l~]f- ^2n, 



n=0 



(2n)! 



and ||\&Q 11^ = |(ioP S^o (tan(6'))^" ' (2nY\" • '^^ conclude that Assumptions 1 and 2 are verified, 
we still have to verify that the vectors we have found here belong to D°°{Bg) and D°°(Al) 
respectively. This would produce a rather long computation, which can be made much simpler 
if we work directly in the coordinate representation, i.e. considering the expressions (14. 3 p of 
our pseudo-bosonic operators. From now on, this will be our point of view. 
Equation Aocp^ = becomes (e*^a; + e~*^^) ^Pq (x) = whose solution is 



cp\,">{x) = N,exp{--e'^'^x'}, (4.4) 



where A'^i is a normalization constant. Analogously, Bq^^q = becomes (e '^x + e*^^) \I'g 



t,T/W - n Koo^v^oo (^-ie^ I ^ied\ .tM. 



X) 



0, so that 



^f\x) = iVsexp |-ie-2'^a;4 , (4.5) 

where, again, N2 is a normalization constant. Notice that, since 9fJ(e^^*^) = cos(26') > for all 
9 E I, both LfQ (x) and \I^q (x) belong to £^(]R), which is the Hilbert space 7i of the theory. 

Defining now the vectors ipn (x) and \l/n (x) as in (12.21) . we find the following interesting 
result: 



V'^\x) = ^ B^ ^f{x) = ^H^ (e^^x) exp {-i e'^' x'} , 
'^^n\x) = ^ {Alr¥^\x) = ^H^ (e-^'x) expj-ie-^^^x^: 



(4.6) 



where Hn,{x) is the n-th Hermite polynomiaL These equahties can be proved using induction 
on n and the following identity of Hermite polynomials: Hn+i{x) = 2xHn{x) —2nHn-i{x). The 
norm of the vectors in (I4.6p can be computed using the following formula: 

2"~^ n! a/vt ,,r, „ .„/r, „ I be 



/°° 9 2"~-^r)' /tt 

e-P^ Hnibx) H4cx) dx = ^^^^;j^ ib' + c'- pY'^ P, 



a/j»(6^ + c^ — p) 



which holds for all p with positive real part, [12]. Here P„ is the n-th Legendre polynomial. 
Therefore we find 



and 



W||2 = |iVJ2 cos (^^ ^ P ( i ^ 



i^i^)ii^ = iiV2rcos(^^]p/ ^ 



^cos(2^)y Vcos(2e) 

which are both well defined (even if the argument of P„ does not belong to the interval [—1,1]), 
for all fixed n. Hence Assumptions 1 and 2 are necessarily satisfied. This was already clear from 
(14. 6p . since the product of a polynomial of any order times a gaussian is square integrable. What 
is not clear at this stage is whether the sets J-"^ = {</)„ (x), n > 0} and J-"^ = {\E'n (x), n > 0} 
are (i) complete in £^(ffi); (ii) Riesz bases. 

To answer to the first question we use the same general idea adopted in [3], which was 
based on a result discussed in [13]: if p{x) is a Lebesgue-measurable function which is different 
from zero almost everywhere (a.e.) in M and if there exist two positive constants 5, C such that 
\p{x)\ < Ce~^l^'l a.e. in M, then the set {x"p(x)} is complete in £^(]R). 

Having this in mind, first of all we notice that J-"^ is complete in £^(]R) if and only if 
the set J^n := {iTn (x) = x" exp { — | e^*^ x^ } , ra > 0} is complete in £^(M). Hence, because 
of the above cited result and since exp | — | e^*^x^} satisfies for our values of 9 the conditions 
required to p(x), then Tn is complete and, as a consequence, J-",^ is complete in £^(M). The 
same conclusion can be deduced for the set J-'^, , which is therefore also complete in £^(]R). 
Therefore, Assumption 3 is satisfied: Ti^ = Ti^ = Ti. 

Let us now go back to the biorthogonality of the two sets J^^ and J^\, . Condition 
(y^Q , \E'q \ = 1 is ensured by requiring that A''! A'2 = ^-j^. Hence, with this choice, we know 



that (</?n , '^m ) = 5n,m which can be written explicitly as 

f Hn {e-'^x) H„, (e-'^x) e"^"'""' dx = 5n,m V2"+-7rn!m!, 
which is a non trivial integral which can also be found in [12] . 
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To understand whether our biorthogonal sets are Riesz bases or not we will proceed as in 
the previous section, by introducing a certain unbounded operator and showing for instance 
that this operator maps an o.n. basis of T-L into the set J-",^ . This would imply that J-l^ cannot 
be a Riesz basis. 

Let us introduce the following unbounded, self-adjoint and invertible operator Tg = e^^^"" ~" •'. 
Then we have 

Ae = TeaTs\ Be = Tea^Ti\ (4.7) 

This implies that Hg = TghgT^^^ where hg = tug (^a'^a + |l). Hence, as in the previous model, 
Tg is an 10: 

HgTg = Tehg, TgHl = hgTg (4.8) 

The same arguments discussed in Section III show that, if the eigenvalues Un = uj{n + 1/2) 
are non degenerate, then a single complex constant a must exist such that 

^l^^=aTg^r^, and ^(f) = =T,-Vn (4.9) 

a 

These equalities show, in particular, that neither J-"^ nor J-"^ are Riesz bases. Also, we 
deduce that S",}, = |apT| and S\, = \a\~'^TQ'^. This is in agreement with the following 
(formal) computations: 

oo / oo \ 

n=0 \n=0 / 

as well as 



n=0 \n=0 / ^ ^ ' ' 



_Tn-2 _ q{0) 



The resolution of the identity looks like 



oo / oo \ ^ 

E ki'^)(^f I = «^^ E \^r.)M ^T- = 1. 
n=0 \n=0 / 



Before closing the section it is interesting to notice that these results could be slightly 
generalized by reversing the point of view we have considered so far: up to now we have 
considered non self-adjoint hamiltonians which can be written essentially as pseudo-bosonic 
number operators. Then we have discovered that an unbounded map exists which transforms 
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the pseudo-bosonic operators into standard bosons, see f l3.8p and (14. 7p . Now we start considering 
the unbounded map Ta^p = e"" "''^" depending on two real parameters a and (3 and we 
define two operators Aa^fs = Ta^fjaT'^^ and Ba^js = Ta^isa'^ T~^^. It is possible to deduce that 

A.a,i3 = a cos(-\/4a/3) — \/f sin(A/4a/3) a"!" and B^^fs = a'^ cos(-\/4a/3) + . /| sin(A/4a/3) a, which 

is clearly different from A^^o, in general. Moreover, [Aa^fs, B^./s] = 1. Hence, we have easily 
constructed a pseudo-bosonic commutation rule. However, there is not a big difference between 
these operators and the operators Aq, Bg considered above in this section. For this reason we 
will not construct the biorthogonal sets arising from A^^is and B^^/s'- the main steps will not 
differ significantly from those already considered here. 

V Conclusions 

In this note we have shown that some non self-adjoint hamiltonians related to relevant quantum 
models can be analyzed within the framework of pseudo-bosons. It is also shown that these 
models give not rise to regular pseudo-bosons, since the biorthogonal bases they produce are 
not Riesz bases. This result, together with other models considered so far, suggests that 
pseudo-bosons may have a real physical interpretation while regular pseudo-bosons are mainly 
mathematical objects. This suggests that unbounded operators play a crucial role in the analysis 
of physical pseudo-bosons, which, in a sense, is in agreement with all the literature on the 
subject: most quantum mechanical systems (but for those living in finite dimensional Hilbert 
spaces, for which the operators are just matrices whose spectra are automatically bounded) 
depend on unbounded operators! Regularity looks more like a mathematical (rather than 
physical) requirement. Consequently, the role of Riesz bases appear more interesting from a 
mathematical rather than from a physical point of view. 

A natural question to consider is now the following: how much of this general structure can 
be extended to generic non-hermitian hamiltonians? This is work in progress. 
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